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Abstract 

Dissipations, e.g. heat flow and bulk and shear viscosities, cause the transport of energy, 
momentum and angular momentum, which is the essence of accretion of matters onto celestial 
objects. Dissipations are usually described by the Fourier and Navier-Stokes laws (classic 
laws of dissipations). However the classic laws result in an infinitely fast propagation of 
dissipations. In relativistic formulation, the classic laws of dissipations violate the causality, 
and hence no relativistic theory of accretion flow onto celestial object is formulated. In 
this short report, we summarize the causal dissipative hydrodynamics, so-called Extended 
Irreversible Thermodynamics (EIT), with a supplemental comment of which the original 
works of EIT are not aware, and then show two theorems about relativistic dissipative flows 
around a Schwarzschild black hole. By these theorems, a significant property of EIT in 
contrast with classic laws of dissipations is clarified, and a dissipative instability of an exact 
solution of relativistic perfect fluid flow is also obtained. 

1 Extended Irreversible Thermodynamics (EIT) 

We begin with summarizing the basis of perfect fluid and classic laws of dissipations (e.g. Fourier and 
Navier-Stokes laws) in order to clarify the basis and need for EIT. The perfect fluid is a phenomenology 
assuming the local equilibrium of fluid, which requires that each fluid element is in a thermal equilibrium 
state. By the local equilibrium assumption, any fluid element in perfect fluid evolves adiabatically and 
no entropy production arises in the fluid clement. This contradicts the dissipative phenomena which 
are essentially the irreversible and entropy producing processes. Therefore, the basic equations of perfect 
fluid can not include any dissipation. Then we notice that the local equilibrium assumption is inconsistent 
with the irreversible nature of dissipative phenomena. In other words, dissipations can not exist in 
local equilibrium systems. On the other hand, recall that the classic laws of dissipations are also the 
phenomenologies assuming the local equilibrium. Hence the classic laws lead inevitably an unphysical 
conclusion; an infinitely fast propagation of dissipations. While the infinitely fast propagation may be 
harmless to Newtonian theories, however, in relativistic theories, it gives rise to a serious problem; the 
violation of causality. The classic laws of dissipations can not be accepted as basic laws of relativistic 
dissipations. (See [1] for details of the violation of causality by the local equilibrium assumption.) 

From the above, it is recognized that we should abandon the local equilibrium assumption in order to 
obtain a consistent dissipative hydrodynamics. Therefore the idea of local non- equilibrium is necessary, 
which means that the fluid element is in a non-equilibrium state. A physically consistent phenomenology 
of dissipative hydrodynamics, which is based on the local non-equilibrium idea, is already formulated. It 
is called the Extended Irreversible Thermodynamics (EIT). As precisely explained in pQ, because of the 
local non-equilibrium idea, EIT describes the entropy production inside each fluid element, which results 
in a finite speed of propagation of dissipations in both non-relativistic and relativistic situations. The 
relativistic EIT is a causally consistent dissipative hydrodynamics [2j El SI IS] - 

Although the EIT is a dissipative "hydrodynamics" , it is called "thermodynamics" . This name puts 
emphasis on the replacement of local equilibrium idea with local non-equilibrium one, which is a rev- 
olution in thermodynamic treatment of fluid element. The terminology "EIT" is found in [1] which is 
developed by experts in non-equilibrium physics. Contrary, the original works of relativistic dissipative 
hydrodynamics [2 , 3 , 4] [5] put emphasis not on the thermodynamic revolution but on the preservation 
of causality. Here let us dare to use the term "EIT" since the local non-equilibrium nature of dissipative 
fluid is the physical origin of preservation of causality. 
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Before showing the basic equations of EIT, we list the basic quantities; 



^(x) : velocity field of fluid (four- velocity of fluid element) 

A mass density (p) and specific volume (volume per unit mass, V) 



V(x) 

Pnc{x) 

T nc {x) 
q»(x) 
U(x) 

n^{x) 



non-equilibrium specific internal energy (internal energy per unit mass) 

non-equilibrium pressure 

non-equilibrium temperature 

heat flux vector 

bulk viscosity 

shear viscosity tensor 

spacetime metric 



Here x denotes the dependence on spacetime point. All of the above quantities except for u M and g M „ are 
the non- equilibrium thermodynamic state variables which characterize the non-equilibrium state of each 
fluid element. The non-equilibrium state variables are classified into two categories: The quantities e nc , 
Pnc, T ne and p are the state variables which exit even at the local equilibrium limit, while the quantities 
<7 M , II and II ^ are the state variables which should vanish at the local equilibrium limit. We call the 
first category (e nc , p n c, T nc , p) the non- equilibrium scalars, and call the second category (q^, II, H 11 ") 
the dissipative fluxes. The dissipative fluxes are the origin of dissipative phenomena and make the fluid 
being non-equilibrium. The suffix "ne" for non-equilibrium scalars denotes "non-equilibrium" , and the 
variables with this suffix have different value from an equilibrium case. Note that the definition of p 
(or V) is the same for both equilibrium and non-equilibrium cases; to count the number of composite 
particles or measure the mass per unit volume. So the suffix "ne" is not given to p and V. 

As the basic assumption of EIT, it is required that the dissipative fluxes are independent non- 
equilibrium state variables. Furthermore, concerning the non-equilibrium scalars, it is also assumed 
that, as in the ordinary equilibrium thermodynamics, the number of independent non-equilibrium scalars 
is two for closed systems which conserve the number of composite particles, and three for open systems 
in which the number of composite particles changes. The remaining state variables are expressed as 
functions of the independent variables through the equations of state, e.g. the non-equilibrium specific 
entropy s no is a function of independent state variables s nc = s n c(£ne, V, g M , II, II MI/ ), where e nc and V 
are chosen as independent non-equilibrium scalars. 

Here note the fact that EIT, at least in its present status, is not necessarily applicable to any non- 
equilibrium state of dissipative fluid pQ. The consistent basic equations of EIT can be formulated for 
sufficiently weak dissipative fluxes, which means that the strength of non-equilibrium nature should not be 
so strong. For example, it is shown in [5] that the EIT preserves causality of heat flow if ^Jq^q^j 'e no < 0.08 
for non-viscous fluid (n = 0, II'"' = 0) under stationary and homogeneous (spatially one-dimensional) 
condition. We can recognize that the amount of energy transported by dissipative fluxes should be less 
than a few percent of the internal energy which includes mass energy. This efficiency of dissipative energy 
transfer is larger than the efficiency of hydrogen burning in a star (~ O(10 -3 ) ). Therefore, although 
the EIT is applicable to a dissipative fluid whose local non-equilibrium states are not so far from local 
equilibrium states, it is expected that the EIT is applicable to many astrophysical systems. 

In the present status of EIT, the above restriction is reflected in the equations of state. The equations 
of state, e.g. s ne (e n c, V, q^ 1 , n, n Ml/ ), are expanded up to second order of dissipative fluxes about the 
fiducial equilibrium state. Here the fiducial equilibrium state is defined as an equilibrium state of fluid 
element of imaginary perfect fluid possessing the same fluid velocity u^(x) and mass density p{x) with 
our actual dissipative fluid. Then, as explained in [T], the evolution equations of dissipative fluxes are 
obtained under two requirements; (i) positivity of entropy production rate at each fluid element, and 
(ii) consistency with experimentally determined phenomenology of relaxation processes of dissipations. 
The resultant evolution equations are; 

r h q» = -q"- \Tu» + T h (q»u„)u»~\A^ ( - T 2 ( /3 hb n , + /3 hs ft )) (1) 

r b ri = -n-c^ + AbCiV^ (2) 

T s (n^)' = -U^ +2T s u a IL ai - v u u) -2r,\u^ v -T (3 hs q^ v ]\ (3) 
where T is the temperature of fiducial equilibrium state. Here, definitions of mathematical symbols are; 
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Q := vfQw, A"" := u^u" + g^ and [A^]° := A^A^A^ - {l/3)A» v A af3 A aP . And the meanings of 
phenomenological coefficients are; r's are relaxation times of dissipative fluxes, A is heat conductivity, £ 
and 77 are respectively bulk and shear viscous rates, and /?hb and /3hs are interaction coefficients between 
<7 M and II or H^. For example, (3^ means that the heat flow arises in a shear viscous flow, also shear 
viscosity arises in a flow with heating. Also /?hb means the same between g M and II. The values of 
those phenomenological coefficients should be determined by some experiments. Furthermore note that, 
because dissipative phenomena are not time reversal, these evolution equations of dissipative fluxes are 
also not time reversal. 

In reference [3] which is one of the original works of relativistic EIT, the molecular viscosity are 
assumed as the physical origin of viscosities II and II ^ v . However in Newtonian accretion disk theories [6] , 
the angular momentum in a disk is transported by the turbulent viscosity. It seems to be appropriate to 
keep various possibility of the physical origin of viscosities. Hence, we treat the viscosities n and n w as 
phenomenological variables, whose physical origin is not specified. This means that the bulk and shear 
viscous rates, £ and 77, are left as undetermined parameter in the equations (|XJ) ^ (ED - As mentioned 
in previous paragraph, £ and rj are the phenomenological parameter whose values should be determined 
empirically. 

The others of EIT's basic equations are given by the conservation of mass current, (pu^) ;At = 0, and 
that of stress-energy-momentum tensor, T^ h = ; 

p + pu».„ = (4) 
p(s+pV) = -q^-q^u^- (nA^+n"") (5) 

(pe + p + Il) ii» = -q» + q a u a u»-u a . a q»-q a u». a -A» a ( (p + II), Q + fi Q % ) , (6) 

where e, p and T are the state variables of fiducial equilibrium state. The above equations (JTJ) ~ (JSJ) and 
the Einstein equation are the basic equations of EIT. 

Finally let us make a comment of which the original works of EIT (TJ [2j [3J [4l [5] are not aware: The 
EIT can not be applied to radiation fluids as shown in [7] . Non-equilibrium radiation fluids need a special 
treatment different from the other matters. However, although the inconsistency of EIT with radiation 
fluid has been revealed, no satisfactory non-equilibrium phenomenology of radiation fluid exits at present. 



2 Dissipative Accretion Flow onto a Schwarzschild Black Hole 

As a preliminary report, we show two theorems, without precise proof, on dissipative flows around a 
Schwarzschild black hole. Both theorems can be proven by solving the basic equations of EIT with- 
out solving the Einstein equation and fixing the metric on Schwarzschild spacetime. The use of fixed 
background metric means that our analysis does not include the self-gravity of the dissipative fluid. 
Therefore the theorems shown below are applicable to non-self-gravitating dissipative flows, i.e. to the 
dissipative fluid whose total energy is sufficiently less than the mass of central black hole, however "non- 
self-gravitating" dissipative fluid is usually assumed in Newtonian accretion disk theories [6]. 

For the first, let us consider the spherically symmetric dissipative accretion flow onto a Schwarzschild 
black hole: 

Theorem 1 Bulk and shear viscosity do not vanish in stationary spherically symmetric dissipative ac- 
cretion flows onto Schwarzschild black hole. 

Naively one may think that the shear viscosity should vanish in spherical accretion flow. However 
theorem 1 denies this naive sense. Here note that the Navier-Stokes law also allows the existence of 
shear viscosity in spherical flow [8]. The difference between EIT and Navier-Stokes is that, while a flow 
wit only bulk (or shear) viscosity is possible in Navier-Stokes law, but such flow in EIT has an unphysical 
divergence; the density or speed of accretion diverges in the region far from black hole. Hence, in EIT, the 
bulk and shear viscosities do not arise separately, and they interact always. This result about viscosity 
may be a significant property of EIT in contrast with Navier-Stokes law. Now I am trying to extend 
this theorem to a non-stationary spherical dissipative flow on general non-stationary spherical spacetime 
which includes the self-gravity of dissipative fluid. 

For the second, let us consider the stationary toroidal dissipative flow around a Schwarzschild black 
hole: 
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Theorem 2 For stationary toroidal flow with no mass accretion onto black hole, any flow of finite 
temperature is impossible, but only a rigid toroidal rotation of zero-temperature is possible. The fluid 
velocity and dissipative fluxes of the possible zero-temperature flow are 

u* = , u r = u e =0 , u' p = nu t and g" = II = 5"" = 0, (7) 

y7(r) - Q 2 r 2 sin 2 9 

where (t, r, 9, <p) is the Schwarzschild coordinate, f(r) := 1 — 2M/r with M as the black hole mass, and Q 
is the angular velocity measured by a rest observer. Since fl is constant and dissipations vanish even when 
they are switched on, the velocity w M in equation expresses a rigid rotation. And the thermodynamic 
state variables of this flow are 

T = p = e = , p(r, 9) = arbitrary in the region f(r) > f2 2 r 2 sin 2 9 (8) 

Here recall that the internal energy e(x), for a fluid element at x, includes the mass energy, kinetic energy 
and the potential energy by external gravity due to black hole. Therefore the vanishing internal energy 
e = means that, for each fluid element, the potential energy due to black hole cancels out the kinetic 
and mass energies of fluid element. 

This theorem together with the third law of thermodynamics implies that any flow of finite temperature is 
non-stationary and/or non-toroidal, and never relaxes to the flow given in equations and (j8]). Hence, 
although a differentially rotating and stationary toroidal flow of finite temperature is possible for "perfect" 
fluid [9] , such stationary flow becomes unstable once the dissipations are switched on. Theorem 2 implies 
the dissipative instability of an exact solution of relativistic perfect fluid. 

The dissipative instability, the EIT's effect, may describe a formation of out-flow from an accretion 
disk without introducing a magnetic field (so-called MRI) , while the "acceleration" of out- flow to form a 
jet will be explained with magnetic fields. In (near) future, it is expected to apply the EIT to a plasma 
composing an accretion disk. 

References 

[1] Jou D., Casas- Vazquez J. and Lebon G., 1988, Extended Irreversible Thermodynamics, 
Rep.Prog.Phys.51, 1105. 

Jou D., Casas- Vazquez J., Lebon G., 2001, Extended Irreversible Thermodynamics (3rd Edition), 
Springer Verlag, Berlin. 

[2] Israel W., 1976, Nonstationary Irreversible Thermodynamics: A Causal Relativistic Theory, 
Ami.Phys.100, 310. 

[3] Israel W. and Stewart J.M., 1979, Transient Relativistic Thermodynamics and Kinetic Theory, 
Ann.Phys.118, 341. 

[4] Hiscock W.A. and Lindblom L., 1983, Stability and Causality in Dissipative Relativistic Fluids, 
Ann.Phys.151, 466. 

[5] Hiscock W.A. and Lindblom L., 1988, Nonlinear pathologies in relativistic heat-conducting fluid 
theories, Phys.Lett.Al31, 509. 

[6] Kato S., Fukue J. and Mineshige S., 2008, Black-Hole Accretion Disks, Kyoto Univ. Press. 

[7] Essex C, 1984, Radiation and the violation of bilinearity in the thermodynamics of irreversible pro- 
cesses, Planet. Space. Sci. 32, 1035. 

Essex C, 1990, Radiation and the Continuing Failure of the Bilinear Formalism, Advances in Ther- 
modynamics 3, 435. 

Saida H., 2005, Two-temperature steady state thermodynamics for a radiation field, Physica A356, 
481. 

[8] Ray A.K., 2003, Viscosity in spherically symmetric accretion, Mon.Not.R.Astron.Soc.344, 1085. 

[9] Fishbone L.G. and Moncrief V., 1976, Relativistic fluid disks in orbit around Kerr black hole, Astro- 
phys.J.207, 962. 



4 



